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A  priori  estimates  for  solutions  of  the  quasilinear  hyperbolic-parabol ic 
equations  governing  the  initial  value  problem  describing  the  motion  of 
compressible,  viscous  and  heat-conductive,  Newtonian  fluids  are  derived  by 
means  of  a  new  energy  method.  This  technique  enables  us  to  simplify  and  unify 
our  previous  results  on  the  global  existence  in  time  and  uniqueness  of  smooth 
solutions  of  these  equations  for  sufficiently  smooth  and  "small"  initial  data 
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SIGNIFICANCE  AND  EXPLANATION 


The  motion  of  compressible,  viscous  and  heat-conductive  Newtonian  fluids 
is  described  by  a  system  of  partial  differential  equations  which  is  of 
hyperbolic-parabolic  type  and  highly  nonlinear.  One  of  the  first  mathematical 
problems  associated  with  this  system  is  the  initial  value  problem.  We  obtain 
the  existence  of  a  unique  smooth  global  solution  in  time  for  the  initial  value 
problem  and  also  the  decay  rate  of  the  solution  as  time  tends  to  infinity. 
Since  the  system  is  quasilinear  with  respect  to  the  unknowns;  density, 
velocity  and  temperature,  we  need  to  assume  that  the  initial  data  are  close  to 
the  constant  equilibrium  state.  The  purpose  of  this  paper  is  to  obtain  a 
priori  estimates  for  the  solutions  of  these  equations  by  means  of  a  new  enerqy 
method.  This  technique,  although  still  necessarily  laborious,  enables  us  to 
simplify  and  unify  our  previous  results,  described  briefly  in  the  abstract  and 
obtained  jointly  with  T.  Nishida  (see,  e.g.  MRC  TSR  #1991). 


The  responsibility  for  the  wordinq  and  views  expressed  In  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


AN  ENERGY  METHOD  FOR  THE  EQUATIONS  OF  MOTION 
OF  COMPRESSIBLE  VISCOUS  AND  HEAT -CONDUCTIVE  FLUIDS 


Akitaka  Matsumura 

$1.  Introduction  and  Main  Theorem. 

In  previous  papers  [1]#  [2] ,  we  have  investigated  the  global  solution  in 

time  of  the  initial  value  problem  for  the  following  equations  governing  the 
motion  of  isotropic  Newtonian  fluids; 


(1.1) 


pt  +  (pU  ’x.  "  0 


U1  +  +  —  p  =  — fy(u^  +  )  +  vi'u  5^)  *  i 

t  x  o  x.  0^  X.  X.  x  ;X. 

j  !  D  1  K  3 


0P 


1,2,3, 


+  u^9  +  -  uj  =  -  (  k9  )  +  40  , 

x .  pc  x .  pc  x .  x . 

3  H  V  3  3  3 


with  the  initial  data 


(1.2) 


(p,u,Q)(0,x)  =  <  PQ 'u0' %  >  • 


3  12  3 

where  t  >  0,  x  =  ^X1'X2,X3^  e  R  •  P  is  the  density,  u  =  (u  ,u  ,u  )  is  the 

velocity,  9  is  the  absolute  temperature,  p  »  p(p,9)  is  the  pressure, 

U  “  u(p,0)  an d  p'(p,9)  are  viscous  coefficients,  <  =  x(p,9)  is  the 

coefficient  of  heat  conduction,  c^  =  cv(p,0)  is  the  heat  capacity  at 

constant  volume,  and  f  =  ~  +  uk  )2  +  \i *  (u^  )2  is  dissipation 

2  x.  x ,  x , 

K  3  3 
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function.  We  shall  consider  the  solutions  only  in  a  neiqhbourhood  of  any 
fixed  constant  state  (p,u,9)  *  (p,0,G)  where  p,0  are  any  positive 
constants.  Moreover,  we  shall  make  the  following  natural  assumptions  on  the 
hyperbolic-parabolic  system  (1.1)  throughout  this  paper? 

(i)  p,  c^,  \i,  p'  and  <  are  smooth  functions  of  0  »  {(p,u,0)  :  |p  -  p| , 

I u | ,  | q  -  e|  <  yq)  / 

2 

(ii)  p  ,  p  ,  c  ,  u,  <  >  0  and  p'  +  r  u  >  0  in  0, 
p  0  V  3 

where  <  min(p,9). 

In  [1],  we  succeeded  in  obtaining  a  global  solution  in  time  of  (1.1), 
(1.2)  by  using  energy  methods  which  were  rather  technical  and  complicated 
under  the  assumptions  that  the  fluid  is  an  ideal  and  polytropic  gas  and  that 

(p0  "  P'  Uq,  0q  -  6)  is  suitably  small  in  H3  , 

where  H  reoresents  the  usual  Sobolev's  spaces  with  the  norm  II  •  R  .  We  also 
proved  the  decay  of  the  solution  to  the  constant  state  ( p,  0,  9),  but  we 
were  not  able  to  estimate  the  decay  rate  of  the  solution.  In  [2] ,  we  obtained 
both  the  global  solution  in  time  of  the  original  problem  (1.1),  (1.2)  and  its 
decay  rate  to  the  constant  state  by  using  a  calculation  of  the  decay  rate  for 
the  linearized  equations,  together  with  the  energy  estimates.  Moreover,  we 
had  to  investigate  the  precise  properties  of  the  spectrum  of  linearized 
equations,  and  assume  that 

—  —  *ii 

(pg  -  p,  u  ,  0^  “  0)  is  suitably  small  in  H  n  L 


The  purpose  of  this  paper  is  to  employ  a  different  energy  method  to 


handle  the  nonlinearity.  The  present  approach  is  simpler  and  less  technical 
than  the  energy  form  used  in  [1],  moreover  we  do  not  require  all  of  the 
precise  properties  of  the  linearized  equations  obtained  laboriously  in  [2] . 
However,  we  note  that  because  of  the  relative  roughness  of  the  new  method,  the 
coefficients  of  the  various  estimates  might  be  more  rough  than  those  in  [1] 
and  [2].  By  making  use  of  this  approach,  we  obtain  both  the  existence  of 
global  smooth  solution  in  time  as  well  as  its  decay  rate  to  the  constant  state 
for  the  general  case  (1.1),  (1.2)  under  the  assumption  that 


—  —  3 

( -  p,  Ug ,  9q  -  9)  is  suitably  small  in  H  , 


(the  previous  approach  required  H4  n  L1).  This  method  will  be  applied  to  an 
initial  boundary  value  problem  for  (1.1),  (1.2)  in  a  subsequent  paper. 

The  main  result  is: 

Theorem  1.1.  Consider  the  initial  value  problem  (1.1),  (1.2)  and  suppose 
the  initial  data  ( pQ  -  p,  u0,  0Q  -  0)  e  H3.  Then  there  exist  positive 
constants  and  CQ  such  that  if  Wp^  -  p,  u^,  0^  -  the  problem 

(1.1),  (1.2)  has  the  unique  global  solution  in  time  satisfying 

p  -  p  e  C°(0,  +«;  H3)  n  0^(0,  +«•  H2)  , 

(u,  0  -  9)  e  C°(0,  +»;  H3)  n  c'to,  +oo*  H1)  , 


sup  I ( p-p,  u,  0  -  0) ( t) 


<  cQn+t) 


3 

4 


lpn-P'  u, 


0  9 


°o  - 


on. 


(1.3) 


X 


$2.  Proof  of  Theorem 

First  rewrite  the  problem  (1.1),  (1.2)  by  the  chanqes  of  variables: 

P  +  p  +  p,  u  +  u,  0  +  0+9,  p( p+p,  0+9)  +  p(p,0),  p(p+p,  9+0)  ->  u(p,0)  and 

so  on  resulting  in 


(2.1) 


(2.2) 


Pt  +  (p  +  p)u3  +  ujpx  -  0  , 


u+uu  +pp  +  -  9 

t  x .  p  x .  —  x 

3  i  P+P  i 


—  (pfu1  +  uj  )  +  U'uk  6i;)) 

—  y  x ,  X.  X 

P+P  j  i  K  j 


i  <e+0)P 

e.  +  uJ9  +  — - uJ  - -  ((<e  )  +  y)  , 

j  (P+P)cv  Xj  (p+p)cy  xj  Xj 


( p, u, 9) ( 0,x)  »  (pQ/  u  ,  0  )  ( x ) 


Define  a  positive  constant  Eq  bv  the  Sobolev's  lemma  so  that  for 

Hf  II ^  <  Eq  we  have  sup|f|  <  CUffl^  <  •  Denote 

a1  ol  a. 

nr  *  {9  f/9x^  3x^  9x^  for  all  a,  =  k)  and  define 

II  •  II  a  II  •  II  •  Then  solution  of  (1*1),  (1.2)  is  sought  in  the  set  of 

functions  X ( 0 ,  +®;  E)  for  some  E  <  E^ ,  where  for  0  <  <  t  <  +«, 

X(tr  t2s  E)  -  {(p,u,0)  ;  p  £  C°(t1,  t2;  H3)  n  t2;  H2)  , 

n^j  c  L2(tl(  1 2;  H2),  (1+t)1/2D2p  £  L2  ( 1 1  ,  t  ;  L2 )  , 

( 1+t)D3p  £  L2(t1,  t  ;  L2)  , 

!u,3)  c  C°(tl(  t2i  H3)  n  chtj,  t2;  H1)  , 


2  3 

fi(u,0)  £  L  (t  ,,  t  ;  H  )  ,  (  1+t) 


>V(u, 


2  2 
9)  £  L  (t  ,  t2;  L  ) 


(  1+t )  H 3  ( u  ,  0 )  £  L2(t1  ,  t  ;  H1  )  , 

anrl  N(t,,  t-,)  <  E  (E  <  E  )}. 
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We  define  N(t1#  t2)  by 

N2(t  ,  t  )  *  sup  (  ll(p#u,  0)  (t)  II2  +  t  MD(  p,u,  0)  ( t)  H2  +  t2  • 

.  HD2(p,u,0)(t)  «^)  +  /t2  iid p ( t )  ii 2  +  nD(u,e)(T)n2  + 

+  t iid2 { p#u, 0) ( t)  ii 2  +  t2(»d3p(t)»2  +  HD3 (u, 0) ( t)  B2) dT 

We  prove  Theorem  1 . 1  by  a  combination  of  a  local  existence  result  and  a  priori 
estimate  for  the  solution  in  X. 

Theorem  2,h  (local  existence)  Consider  the  initial  value  problem  (2,1) 
for  t  >  t^  with  the  initial  data  at  t  *  t^  a£ 

(2.3)  (P'U'OHt^  c  H3  . 

Then  there  exist  positive  constants  <  EQ)  and  t  which  are 

independent  of  t^  such  that  if  N(t^,  t^)  <  ,  the  problem  (2.1)-(2.3)  has 

the  unique  solution 

<p,u,G)  e  x(t1#  t1  +  t?  C1N(tl,t1))  . 

The  proof  of  Theorem  2.1  is  given  in  the  same  way  as  in  [1]. 

Theorem  2.2.  (a  priori  estimates)  Suppose  that  the  initial  value 
problem  (2.1),  (2.2)  has  a  solution 


(p,u,0)  £  X( 0 , T;  E) 


for  some  T  >  0  and  some  E  <  EQ.  Then  there  exist  three  positive  constants 

/  2* 

e2  and  e3  ^e2'e3  *  and  C2  ^C2e3^1+Cl  *  Eo^  w^lch  are  independent  of 

T  such  that  if  **po#uo'®o'*3  *  G2  anc*  E  <  g^,  then  the  solution  satisfies 
the  a  priori  estimate 

(p,u,0)  e  X(0,T;  c2  HP0  'V90  ll3)  • 

Proof  of  Theorem  2.2.  Take 


e3  G1  e3 

G  =  min(£l,G2,  ,  --7""2) 

1  2  C!  /1+C 


We  may  use  the  standard  continuation  arguments  of  a  local  solution  on 
[0,nT],  n  =  1,2,...  to  get  the  global  solution.  In  fact,  by  the  local 
existence  theorem,  the  definition  of  gq  and  the  assumption 
)lp >%,u^ ,  9^  II  <  c^,  we  have  a  local  solution 


:p,u,0)  g  X(0,t;  Cpp  ,V 


5o  "3  > 


By  C  »pn,u  ,9  II  <  C1€q  <  e3  and  11  /Uo '  90  H3  *  C2 '  the  a  priori  estimate 
gives 


(  p ,  u ,  9 )  g  X(  0  ,  t?  II  , u^ 


0rt 
0  3 


Then  bv  II ,  9^  II ^  *  G )  and  the  local  existence  theorem  with 

i,  we  have  again 


(p,u,0)  g  X(t,  2 1 ;  C1C2  1,pQ  ,U0 ' 90  11 3  5 
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I 


Noting  that 


N  ( 0 , 2t )  <  N  (0,t)  +  N  ( t,  2x) 


we  also  have 


(p,u,0)  e  X( 0 , 2tj  C2/1+C1  »PorUo,0oll3) 


/  2  / 2 

Now  by  c2^1+ci  Bpo'U0/90*3  4  C(r1+C1  e0  <  C3/  the  a  Priori  estimate  shows 


that 


(p,u,0)  E  X{0,2ti  C2  ,PO'UO'0O  *3'  * 


Thus  we  can  continue  to  use  the  same  arguments  on  [nx,  (n+Dt]  and 
[0,  (n+1)tl  successively  n  *  2,3,...  .  Finally  the  estimate  (1.3)  follows 
from  Nirenberg's  inequality  [3] 

1  1 

(2.4)  Bup|p,u,e|  <  Clp,u,9n4  KD2  (  p,u,  0)  I!4 

x 


§3.  A  Priori  Estimates 

We  present  here  an  energy  method  to  obtain  a  priori  estimates  for  small 
solutions  of  equations  with  dissipatin.  First  we  rewrite  the  system  (2.1)  so 
that  all  the  nonlinear  terms  appear  at  the  right  hand  side  of  equations; 


\ 


(3.1) 


/  L°  =  p  +  pu^  -  f°  , 
t:  X  j 

i  i  A  ~  i 

L  hu  +  p„ p  +  p„0  *  yu 

t  F1PX,  h2  Xi  x.Xj 


4  i  ^ 

L  5  6  +  p.U  “  K0  *  f 

t  3  X.  x.x. 

3  3  3 


I'  1 14*  1 1 1  ^  ^  sa  f 


-  (  u+u*  )u 


X.X. 

1  3 


where 


P1  =  Po(0,0)/p,  P2  =  p  (0,0)/pf  P3  =  9p0(O,O)/p  Cv(0,0)  , 

V 


p  2  r0' 

\\  *  y(0,0)/p,  y*  =  y’(0,0)/o#  <  *  <(0,0 )/p  0  (0,0)  , 


f°  H  -  pu"5  -  U3p  , 

X.  X. 

3  3 

f*  =  -  uV  +  (p_  -  p  /(p+p)p  +  (P0  -  p  /(p+p))0  + 

x.  1  p  X.  2  0  x 

j  *  1 

+  (Vx  +  Vx.)(V  +  ux.)/<y  +  p)  + 

+  (u'o  +  u’0  )u3  /(p+p)  +  (u/(p+p)  -  g)u*  + 

p  X ,  0  X .  X  .  *  • 


i  3 


3  3 


+  ((u+u')/(p+p>  -  (y+M,))u 


X  X  . 

1  3 


and 

f4  =  -  u30  +  <P  -  (0+0)p6/fp+p)Cv)u^  +  U  o  +  x09x  )9X 

j  3  3  3 

♦  (</(p+p)Cv  ~  <)0x  +  ^Pcv  • 

3  j 

Define,  for  k,  m  *  0, 1,2,3, 


Ak(t)  =  /  ~  nkf°  •  Dkp  +  Dkfi  •  nV  +  —  nV  •  ok0  dx  , 

P  P3 


m  /  m  0  m  p  _ in  i.  to 

B  (t)  z  f  D  f  •  D  p  +  H —  D  f  •  D  p  dx  , 

Xi  Xi  2u+u'  i 


C*ft)  3  \  2 —  oV  •  D4f°  dx 

2U+U’  *i 


/ ( p+  p ) 

3 


8- 


Then  we  have  the  following: 


Lemma  3 « 1 .  There  exist  positive  constants  v  and  C  which  are  independent 
of  t  such  that 

(i)k  »Dk(p,u,0){t) K2  +  v  Jq  HDk+1(u, 6) (t) n2  dT  < 

<  C(  IIDk(po,uo,0o)  n2  +  |Ak(T)|dT),  0  <  k  <  3  , 

(ii)m  HDmp(  t )  II2  +  v  /q  KDmp(  t)  II2  dt  < 

<  C  C  HDmp0  II 2  +  ID1"  1  N 2  +  HD1"  1  u ( t )  II 2  + 

+  IIDm(u ,  0)  (  t)  II2  +  lBm  1  <  xH  +  |cm  \x)|dx),  1  <  m  <  3 

0  2  ft.  P.  1  2 

(iii)  t  IID  (p,u,e)(t)«  +  v  J  i  IID  (u,0)(t)II  dT 

jL  $  K  v 

<  C(/J  T4-1  IIDk(p,u,0)(T)  n2  +  li£Ak(x)|di),  l  =  1,2,  k  =  1,2,3 
(iv)  t ^  II Dmp ( t )  II 2  +  v  t*  lDmp(  t)  l|2dx  < 

<  C(t*'IIDm_1u(t)  II2  +  /p  T&_  1  IIDmp  (  T )  U2  +  T  ^  1  UDm  1  Xi  (  T  )  U 2  + 

+  i*'IIDra(u,  0)  (  T)  l2  +  I  T^C™  1  (  t)  I  +  I  t^Bm  \x)|dx  , 

l  =  1,2,  k  =  2,3  . 
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Lemma  3 . 2 .  There  exist  a  positive  constant  C  which  is  independent  of  t 
such  that 

2 

(3.2)  N2(0,t)  <  C(  lip  ,u0,9q  II2  +  l  (1+t*)(|A£(t)|  +  |b£(t)|  + 

4=0 

+  |CA( T) | )  +  (1+T2)|A3(T)|dT)  . 


Proof  of  Lemma  3.1. 


For  ( i).  and  (iii) 

*  £,k 


we  may  estimate  the  equality 


f  —  nk(L°)  •  Dkp  +  Dk(L1)  •  DkuX  +  —  Dk(L4)  •  RkQ  dxdt 
0  '  T  P, 


ft  Ik. 

=  J0  T  A  ( T ) dt 


which  implies  after  integrating  by  parts 


4  ,  P1  .  k  .  2  1  .  k  .  2  P2  .  k  . 2  t  t 

t  /  —  Id  p|  +  —  Id  u I  +  - —  |d  9|  dx 

2p  P3  t=0 


rt  Z  r  k+1  .2  k  i  ,2  <P2  4,  k+ 1 . . 2  ,  „ 

+  Jq  t  J  u In  u|  +  (u+u*)|d  |  +  — —  t  Id  0|  dxdi 


P1  .  k  ,2  1  ,  k  ,2  P2  . _k „ . 2 . 


Jn  4t  (i-5  )  J  —  Id  pi  +  -  Id  u|  +  Id  e|  dx  + 


2p 


0  V 

+  T  A  (x)dT  , 


where  fi1  represents  Kronecker's  delta.  These  inequalities  prove  (i)^ 
and  (iii)  easily. 

'C  r  k 

Next  for  (ii)^  and  ( iv) a 

m  £,m 


we  may  estimate  the  equality 


ft  _i  r  .  _m-1  +  p  p”1"1^1)  .  Dm_1 


ft  *>  r  I  .  u  n 

Jn  T  /  D  <Lv  >  *  D 


i  2  ii+11' 


p  dxdx 
xi 


rt  A  m-1, 

J0  T  B  ( T)<*T 


which  implies  after  integrating  by  parts 


t*  /  i^i!  +  „v 

2m+u' 


dx 


t*t 

+ 

t-0 


l 

T 


I-2- 

*  -V  -v 

2g+u’ 


,  _ra  ,2  PP2 

'D  Pi  +  ~~ZrZ~‘ 

2p+y' 


m- 1 

D  P. 


p  .  m-1  i  .2,  „ 

- K —  | D  u  |  dxdx  - 

2H+^'  i 


fn 


fcx  ( 1 


Id"oI2 


m  1 
D  U  - 


m 

D  P. 


dx  + 


2y+y* 


+  C  ( x) t  +  B  ( T) t  dT 


These  equalites  imply  (ii)  and  (iv) 

Proof  of  Lemma  3.2.  For  any  positve  constant  e  >  0,  consider  the  form 


(3.3) 


ll  2  3 

l  (i).  +  e  y  (ii)  +  edit).  .  +  e  (iv) 

k=0  R  m=1  "  1,1 


1 1 2 


*  ‘4|ili>2,2  *  e5|iV>2,3  *  '6|111,2,3  • 


By  taking  c  suitably  small  in  (3.3),  we  can  easily  prove  (3.2)  by  (3.3). 


Next  let  us  estimate  the  nonlinear  terms 


Lemma  3 . 3 


_ .  Suppose  (p,u,6)  e  X(0,T;  E)  for  some  E  <  •  Then/  for  an^ 

positive  constant  c,  there  exists  a  positive  constant  C(e)  which  Is 


independent  of  t  such  that 


2 

/q  l  n  +  |AA(t)|  +  I  B&(  T)  )  +  I C  T)  I  )  +  (  1  +  T2)A3(T)dT 


<  e  N2(0,t)  +  EC(E)N2(0,t) 


Before  proving  Lemma  3.3,  we  note  that  Lemmas  3.2  and  3.3  easily  imply 
the  desired  a  priori  estimates.  In  fact,  we  may  first  choose  e  so  small  and 
next  choose  E  so  small  that  we  have 


N2(0,t)  <  CBp0,u0,Q0H2  for  E  <  e2 


Proof  of  Lemma  3.3.  Because  there  are  many  terms  to  estimate,  we  nick  up  some 
examples.  The  remaining  terms  will  be  estimated  in  the  same  way.  First  let 
us  pick  up 


(3.4) 


Jo  J  f0p 


in  A°(x)dT.  By  using  Nirenberg's  inequality  (2.4)  efficiently  we  estimate 


(3.4)  as  follows; 


f*  f  f°p  dxdil 


I  Jn  J  -  pV 


u  pp  dxdT  I 

j 


=  /  2~  Ul.  dXdTl 

1 
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<  J  /q ( sup |  p | )  II p I  HDutldi 


<  c  Jg  npn4  hduh  «D2pi4dT 


<  CE(sup)pI2  )8(sup(  1  +  t)  HDu«2)2(sup(  1  +  t)2|ID2p|2)8 


£  4  dx 


<  CEN  (0,t)  . 


Next  let  us  consider  the  quantity 


(3.5) 


/q  /  T2  °3(PU^>X  *  °3p  dxdT 


rlnq  in  f*  t2A3(r)dr  or  x2B2(x)dT 


We  estimate  (3.5)  as  follows; 


\ 


rfc  r  VJ  ^ 

Cj  ”j  ‘VjOh  "£"k'm 


(3.5)  *  /  t  (p  +  puD  )  •  p  dxdi 

;0  ;  Kx.  K  x  x.  xx  xx.  x 


rt  r  2  j  3  3  2  j 

J^JluDp  *  D  p  +  3t  p  u  p 

;0  J  x.  K  x.x,  xx  •  x  x  x 

j  3  k  £  m  k  i  m 


2  j  2  3  j  3 

+  3xp  u  p  +TpDu  »Dp  + 

x  X,  XX  X  .  H 

3  k  £  m  k  £  m  3 


2  3  3  1  2  1 

+  t(Dp*Dp)uJ  +  3t  p  uJ 


.  X.  X  ,  X  X  X.  X  > 

j  k  j  t  m  k  2, 


2  £, 

+  3t  p  u  p 

X.  X  .  X  X  X,  X  .  X 

k  3  i  m  k  j  m 


+  X2p  D3u^  •  D3 


x  . 
3 


0  dxdt 
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-  \f0  !  -  2  dxdTl 

<  CE  fj  t2  ID3pU)*2^  <  CEW2(°,t)  , 


.  U7l</5/e,Vol2+^lnV  lo2u,2,xaT 


<  eN2<o,t>  +  \  f^2  (|lD2ol4dK]2(/  |n2u|4dx)2dt 


N2<o,t>  ♦  c  /„  t2«d20«2  «d2u«2  dT 


<  t  N2(0,t)  +  £  (aup( 1+t2 ) *02p#2 )(aup( 1+t2 ) #02ul2 ) 


X  ?»  T20  +  T2)"2dT 


<  e  N2(n,t)  +  Cc'^ee  N2(0,t)  , 


,  |TJ  and  Uft!  are 


estimated  in  the  same  way  as  II-, 


n„l  -  If*  r  r2  P  ■  o’o^  ax«,l 


<  }fc  T2  suplpl  »d3p«  «D4«»dT 

X  3 

<  ft  c  t2„d4uii2  +  C  e' 1 T2  n p#2»d2p»2^3p»2^ 


c  r  M2(0,t)  ♦  CEFa  €~V<0,t>  T2n  +  T> 
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